function Geomatry1D(xmin, xmax, p, K, x)
    """
    Generate geomatry informations.
    """
    h = (xmax - xmin)/K
    VX = [xmin + i * h for i = 0:K]
    @inbounds x = ones(p+1,1)*VX[1:K]'+(x .+ 1)./2*(VX[2:K+1]-VX[1:K])'
    
    return h, x
end